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Abstract 

•i-H . 

■ It is widely believed that at small x, the BFKL resummed gluon splitting function 

should grow as a power of 1/x. But in several recent calculations it has been found to 
decrease for moderately small- a; before eventually rising. We show that this 'dip' structure 
is a rigorous feature of the P gg splitting function for sufficiently small a s , the minimum 
occurring formally at log 1/x ~ 1/ ^foT s . We calculate the properties of the dip, including 
corrections of relative order a/oI, and discuss how this expansion in powers of \foT s , which 
is poorly convergent, can be qualitatively matched to the fully resummed result of a recent 
calculation, for realistic values of a s . Finally, we note that the dip position, as a function 
of a s , provides a lower bound in x below which the NNLO fixed-order expansion of the 
splitting function breaks down and the resummation of small- x terms is mandatory. 



1 Introduction 



A major effort is currently under way to push the precision of DGLAP [1-3] splitting func- 
tions to next-next-to-leading order (NNLO) accuracy [4]. One of the main applications of 
such an effort could be to improve the description of the small- x parton distributions, which 
with the current NLO evolution suffer from pathologies such as negative gluon distributions 
and predictions of a negative Fl [5,6]. Furthermore a good knowledge of small- x parton 
distributions will be ever-more relevant as collider energies are increased, for example at the 
LHC or a possible VLHC, which will be able to probe small- x kinematic regions unexplored 
even at HERA. 

However a question that remains to be understood is that of the domain in which fixed 
order expansions are sufficiently convergent as to be reliable. Indeed it is known that at small 
x, there are large logarithmic enhancements of the splitting function at all orders [7,8], leading 
formally to the breakdown of the convergence of the series for a s log 1/x ~ 1 and it has been 
argued [9] that there is evidence in the data [10] for the presence of some such terms. 

Much effort has been devoted in recent years to resumming these logarithmically enhanced 
terms, which are expected to lead to a rise at small x, as a power of x, for the gluon-gluon 
splitting function, xP gg (x). It turns out however that the LLx summation, a™ log™ -1 1/x rises 
much too steeply [11,12] to be compatible with the more gentle rise of the F2 data [10]. On 
the other hand, the inclusion of the NLLx terms a™ log n ~ 2 1/x — extracted from the NLLx 
kernel eigenvalue [13, 14] and based on several Regge-gluon vertices [15] and on the qq cluster 
[16, 17] - - leads at moderately small x to a negative splitting function [18, 19]. Since that 
discovery, there has been investigation of the origin of these problems, and various approaches 
have been proposed to estimate yet higher orders [20-30], the most successful of them being 
based on a simultaneous treatment of small- x and collinear logarithms. 

A surprising observation, common to all these approaches, is that in the phenomenologi- 
cally relevant, moderately small- x region, the splitting function actually decreases, while the 
power-like rise is delayed to somewhat smaller values of x (resummed curve of figure ^ [23] , 
which has been found to be rather close to a splitting function fitted to the F2 data [27,30]). 
The question arises therefore of whether the resulting 'dip' structure is a well-defined property 
of the gluon splitting function, or instead perhaps an artefact of the particular schemes used 
to 'improve' the small- x hierarchy. The purpose of this letter is to show that the dip has a 
simple origin, specifically in the structure of the first few terms of the perturbative series, 
possibly matched to a resummed behaviour at smaller x values. 

More precisely (Sec. |2J), in the formal limit of small a s , the dip is a consequence of an 
interplay between different fixed orders, and one finds that the simple fixed-order hierarchy 
breaks down not for a s log 1/x ~ 1 as widely expected, but rather for a s log 2 1/x ~ 1. The 
result is that the properties of the dip can be described in terms of a series in powers of y / aj. 
For phenomenologically relevant values of a s though, this series in y/a~ s turns out to be very 
poorly convergent. Instead we find that quite simple resummation arguments, presented in 
section |31 still enable us to gain some quantitative understanding of the dip properties. 

2 Low perturbative orders and ^/c^-expansion 

Let us start by recalling the structure of the LLx terms of the xP gg {x) splitting function, 

A nt n-l <X n s log™- 1 - , (n > 1) , (1) 
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Figure 1: The xP gg (x) splitting function. The resummed (NLLb) curve corresponds to scheme 
B of [23]. 



where a s = a s N c /ir. A number of the lower order terms in the series are absent, A21 = ^32 
= 0, while 
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Since these and all further terms are positive, the LLx splitting function grows monotonically 
as x decreases. The NLLx terms can be written as 



A n , n _ 2 a n s log n ~ 2 - , 

x 



(n > 2) 



where the first few coefficients are [13, 14] 
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and b = ^| — is the first beta-function coefficient. These coefficients are given in the 
Qo scheme [31] and for renormalisation scale fi = Q. They come from a simple expansion 
of the NLLx kernel eigenvalue, and — notably A31 and A42 — can be traced back to early 
calculations of NLLx gluon vertices [15] and of the qq cluster [16, 17]. They include — in 
particular ^42 — the running coupling effects, which are part of the NLLx corrections. In 
the Ms scheme only the rif parts of ^20 an d ^-31 will differ, while from A42 onward the rif 
independent part will differ as well. Because of the zeroes in the LL coefficients, A31 and ^42 
are independent of the choice of \x. 

The resummation hierarchy as written above in terms of LLx and NLLx terms is intended 
to be applied when a^logl/x is of order 1, while q s « 1 and log 1/x 3> 1. Let us however 
examine an intermediate small- x limit in which log 1/x > 1 but a s log 1/x -C 1 (the precise 
region will be better specified shortly). 
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Because the LLx coefficients A21 and A32 are zero, the lowest order term with log 1/x 
enhancement is the NLLx term A^a^ log 1/x, which is NNLO in the usual DGLAP pertur- 
bative expansion. Since A31 is negative it will lead to an initial decrease of the splitting 
function and at some sufficiently small value of x the NNLO gluon splitting function [32] will 
become negative, as shown in figure ^ where we have included the small- x part of the NNLO 
xP gg (x), A\q6i s + Asi6? s log 1/x (A20 = in the particular scheme used in the figure [23]). 

At N 3 LO, order a^, both LLx and NLLx terms are present. Since we are in the regime of 
log 1/x ^> 1, the LLx log 3 1/x term will clearly dominate over the NLLx log 2 1/x term. 
What is interesting however is the interplay between the negative NLLx a 3 log 1/x term and 
the positive LLx log 3 1/x: 

1 1 

xPgg(x) = const. + A 3 i a 3 log - + A43 a s log 3 — H , (5) 

where the constant term includes Aiqo s and ^boo 2 contributions, and at each order in a s we 
have written only the term with the strongest log 1/x dependence. Since ^43 is positive and 
has stronger logx dependence than the negative A31 term, the splitting function as written 
in © will eventually start rising. The A31 and ^43 terms will be of the same order when 
a s log 2 1/x ~ 1, and the splitting function of eq. (JSJ) will have a minimum at 



log— = V-iT Ll - (6) 
x min V 3^43 a s 

The appearance of this minimum for a s log 2 1/x ~ 1 suggests that it may be of use to examine 
an alternative classification of the series, in which we consider all terms that are of similar 
magnitude when a s log 2 1/x is of order one, 1 

A fc , 2fc _ 5 ^log 2fc - 5 i, (3<k<4). (7) 

One finds that there are only terms with k = 3,4, since lower values of k would be associated 
with negative powers of log 1/x, while higher values of k would be super-leading in the usual 
LLx classification. In other words the two terms, a 3 log 1/x and af log 3 1/x, that we have 
examined so far provide the full leading contribution for a s log 2 1/x ~ 1. 

This is illustrated in figure |5J which shows various possible classifications of logarithmi- 
cally enhanced terms. Rows correspond to a given power of a s ; columns to a given single- 
logarithmic order (LLx, NLLx, and so on); terms on a same downward going diagonal line 
(reading from left to right) all have the same power of logx. 

Terms on upward going diagonal lines in figure [21 are of the same order for a s log 2 1/x ~ 1. 
On any given such diagonal, the number of terms is always finite, due to the fact that the 
natural hierarchy is single logarithmic. The leading terms in this regime, discussed above 
in eq. Q, are highlighted by the upper (upward-going diagonal) ellipse. The lower ellipse 
contains terms suppressed insofar as they have one less power of log 1/x, 

^,2fc-6«"l°g 2fe - 6 i, (3<A:<5). (8) 



1 This is a double-logarithmic classification, however one should bear in mind that the perturbative series 
itself contains at most single logarithms — our study of powers of a s log 2 1/x therefore just represents a 
particular way of reclassifying terms in the single-logarithmic perturbative expansion. 
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Figure 2: Representation of different classifications of logarithmically enhanced terms. Sym- 
bols £ x' indicate terms that are present; '0' indicates terms that could have been present but 
are zero; 'nj' indicates a term whose only non- vanishing part is proportional to nf, a dash 
indicates terms which do not exist by definition. 



Equivalently, since we are interested in the region where logl/x ~ l/y / aj, these terms are 
suppressed by a power of y/al- Adding the terms of the lower ellipse to eq. © one obtains 

xP gg (x) = const. + A 31 a 3 log - + A i3 a* log 3 - + A A2 a* log 2 - + O ( a k s log 2h ~ 7 - ) , (9) 

X X X \ X J 

where we have exploited the fact that A54 = and that the ^30 contribution can be absorbed 
into the constant piece. Solving for the minimum of eq. Q gives 



1 1 / ^3i . A\ 2 A 42 nn s 

Xmin V 3A 43 a s 9^| 3 3^ 43 



A 31 1 A 42 . 1.156 



+ 0(^- s )^-— + 6.M7 + 0{V^), (11) 
V 3A i3 a s 3^43 ^a s 

where the numerical values have been given for rif = 4. We see that the effect of the subleading 
A42 term is to shift logl/x m i n by a (rather large) constant. 

As well as considering the position of the dip, it is interesting to study also its depth, d. 
Substituting log 1/x ~ a s into eq. © one immediately sees that the dip's depth is of order 
6^J 2 ■ Including the subleading terms (lower ellipse of figure EJ) gives the following result 



2A 31 f -3A 31 _ 5/2 I A 31 A A2 



A i3 s 3 ^ 43 



a 3 s +0 (a 7 J 2 ^j ~ -1.237a 5 / 2 -11.15a s s +0 (aj /2 ) . (12) 



The depth has been defined with respect to the x = 1 limit of eq. ©, which includes the usual 
6t s constant term, but also A 2 qoi 2 term and the unknown NNLLx term A 3 oa>g. The full P, 



mi 



splitting function has of course a 1/(1 — x) + divergence so its x = 1 value can not actually be 
used as a reference point for defining the depth. So one may choose to define it alternatively 
with respect to the value of the x — > LO splitting function, Aioa s . This introduces extra 
terms A 2 Qa 2 + A 3 q6? s in the expression, eq. (|12j). for — d. 
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Figure 3: Properties of the dip in the NLLb model of [23] compared to our analytical predic- 
tions. See text for details. 



The dip position and depth, as a function of a s , are shown respectively in figures 0k and 
I3J3. In each case the solid line represents the dip properties as 'measured' from the NLLb 
scheme 2 of [23] , which was shown also in figure ^ The shaded band represents the spread 
of the predictions based on eqs. iP|l-p2j). The upper edge of the bands, labelled 'Quadratic 
solution' corresponds to the use of eq. (|TU|) and its direct substitution into eq. ©; the lower 
edge corresponds to eqs. (fTTj) and (fT2*)) . For small values of a s , there is rather good agreement 
between the expanded forms of our predictions and the dip properties as measured from 
the full resummation: the dip position is within the uncertainty band, typically close to 
the expanded solution; the depth is just outside the uncertainty band (again closer to the 
expanded solution), though this may be because we have measured the depth with respect to 
the Aiqoi s reference level and have not included the resulting additional unknown NNLLx ^30 
contribution to the depth. Instead, including the ^30 as it appears in the NLLb model, lowers 
the band so that it overlaps with the measured depth. Leaving aside these details, for both 
the position and depth of the dip, the scaling with a s is clearly reproduced, providing strong 
evidence that the dip truly is a consequence of the low-order behaviour of the perturbation 
series. 

We note though that the spread of predictions, based on eqs. (fTU|) and ((TTj) . is quite 
significant. This is essentially due to the large value of the A42 coefficient, which means that 
the series in ^fal in eqs. (fTT|l and (|T2*|) is very poorly convergent — the leading and subleading 
corrections are of the same order when a s ~ 0.01-0.02. 

In practice our low-order arguments seem to extend somewhat further, providing a rea- 
sonable description of the dip, within the large uncertainties, up to a s ~ 0.05 — 0.1. However 

2 We note that since the NLLb scheme accounts only partially for the n/ dependence (that associated with 
running of the coupling), the resulting NLLs A nt „-2 coefficients differ slightly from those shown in eq. @, 
with A20, A31 and A42 corresponding to the n/ = results of eq. (|1J (in A42 the n/-part in the b-dependent 
term is retained, hence the coefficient of n/ in A42 is 0.334). The reason for the only partial inclusion of the 
rif dependence is that the NLLb scheme is based on a single-channel, purely gluonic approach, whereas full 
account of n/ dependence would require a two-channel, quark-gluon formulation. 
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beyond this point the prediction fails quite dramatically, with the height of the predicted 
dip minimum becoming for example negative (Aioa s — d < 0), in contradiction with the full 
resummed results. Furthermore there is a clear change in the a s dependence for both the 
measured position and depth of the dip. This suggests that for a s > 0.05 the dip description 
can no-longer be founded on low-order perturbation theory alone. 



3 Resummation and cut-representation argument 

On the other hand we know that when a s is moderate and log 1 jx is sizeable we enter 
the usual regime of resummation of terms (a s logl/x) n [7], together with its subleading 
corrections [13,21]. Though the strict LLx, NLLx hierarchy is ill-behaved, the inclusion of 
renormalisation group effects tends to stabilise this hierarchy (e.g. [23,29]). As a result one 
obtains the usual, expected behaviour of a splitting function that increases as a power of x 
at small x. 

A simple estimate of the x value for which this increase occurs can be obtained in the 
approximation of a frozen coupling using the quadratic expansion of the effective BFKL 
characteristic function 

asXeff (7, oc s ) = u s (a s )(l + D(a s )(j ~ 7m) 2 ) , (13) 

where uj s (a s ) is the value of a s \ at its minimum, 7 = j m , and D(a s ) is related to the second 
derivative of x (see figures 1 and 3 of [23]). 

This leads to the well-known square-root branch-point for the anomalous dimension, 



7 = 7 m + A /^— , (14) 



and to the representation 



xP aa (x)~ / — J— x~ w ~ — ■ ^ , (15) 

" y ' J vr V Dlj s log 3 / 2 1/x V ' 

for the splitting function. 

It is amusing to note that the above estimate shows a dip at 

1 3 

u s (a s )log- = (16) 
x 2 

due to the logarithmic prefactor. Of course the actual cut structure of the anomalous di- 
mension is much more complicated, showing a variety of subleading branch cuts, generally 
at complex oj values [11,18], which are needed to match the small-x representation ()15|) to 
perturbation theory for small a s log 1/x. For this reason the dip structure ()16|) . based on the 
moderate-x behaviour of (|15[) is not always to be taken seriously. 3 

3 For example, as we have mentioned earlier, the fixed-coupling LLx splitting function has no dip at all. 
It is interesting also to note that LL evolution with (subleading) running coupling corrections does have a 
dip [26, 33] — its small-a s properties are different from those of the full NLLx dip, because it is due to 
an interplay between terms a^\og n l/x (1 < n < 3) and so, in the limit of small q s occurs for log 1/x of 
order 1. A related running-coupling LLx dip has been obtained in [28,29], though the different scale of the 
running coupling and the use of the Airy extrapolation mean that it has different formal small-a s properties 
from [26,33]. 
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However, in our resummed calculation, the existence of the dip relies on the negative 
log 1/x-slope of the splitting function which is pretty well represented by the y / oj-expansion, 
as noticed before. Furthermore, for iv s logl/x > 3/2, eq. (|15|) is a reasonable representation 
of the splitting function and in cases — as ours — in which there is a dip, we can take eq. ()16|) 
as an upper bound on its position. 

In the running-coupling case it is to be kept in mind that the cut is actually broken up 
into a series of poles, the leading one being at a position uj c {a s ) which lies somewhat below 
io s (a s ) (see figure 18 of [23]) because of running coupling effects. Nevertheless, as long as x 
is not too small, the inverse Mellin transform (|15|) does not resolve the difference between a 
cut and series of poles. 

Therefore, by joining the y / aJ-expansion with the cut-representation arguments, we are 
led to believe that the perturbative and resummed regions can be matched by the inequality 

log ~ + c 2 < —— , (17) 

where c\ and C2 are provided by eq. (|11|). and we have replaced u> s with uj c < uj s . Since the 
right-hand expression goes as l/a s , this equation provides a transition point in a s , below 
which one should use the perturbative (double-logarithmic) representation described before, 
and above which one should use the full resummed behaviour. 

This is confirmed by the moderate a s region of figure where one sees a clear bend 
in the behaviour of logl/x m i n when 3/2u> c becomes of the same order as the perturbative 
representation, eq. (|lljl. with the measured logl/x m i n remaining consistently below 3/2w c . 



4 Conclusions 

The arguments provided in this letter go some way towards explaining the features of the dip 
for a range of a s values, both in terms of a perturbative series in powers of y / aj for small a s , 
and in terms of a resummed upper bound of the dip position, ~ 3/(2u; c ), for moderate values 
of a s . 

It is the moderate-a s region that remains the least well understood, the matching of the 
small- x increase to the initial decrease being a quite non-trivial problem. For example the 
simple resummed treatment given above is subject to additional running-coupling effects (e.g. 
difference between uj c and lo s ) which may contribute further displacement of the dip and which 
have not been considered here. Nevertheless, the arguments given so far show that the dip 
does exist, as a moderately small- £ phenomenon, under the simple condition that the small- £ 
part of xP gg {x) has initially a negative log 1/x slope, as is the case starting at NNLO. 

An important phenomenological point that remains to be made concerns the validity of 
fixed-order expansions of the splitting functions. From our analysis of the dip properties, it is 
clear that for a s > 0.05 one starts to see a breakdown of the perturbative expansion. Despite 
this fact one notes a remarkable property of figure Q namely that the pure NNLO expansion 
of the splitting function coincides rather well with the resummed result up to the position 
of the dip minimum — considerably beyond the point in x where one would have naively 
expected the DGLAP terms to completely change the behaviour of the splitting function. 
This holds for a wide range of a s . 

We cannot claim to have fully understood this observation, however it does suggest that it 
may in general be safe to use the fixed order, NNLO, P gg (x) splitting function down to x values 
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corresponding to the dip position, and only beyond this point will small- x resummation be 
strictly necessary. Thus one can use the 'measured' dip position, the solid curve of figure^, as 
an estimate of the limit of validity of the NNLO expansion at small x. Considering this in the 
context of the available F% data, one sees that the limit cuts through the HERA kinematical 
range, suggesting that while much of the data will be in the region that is 'safe' for an NNLO 
analysis, there is also a substantial region at lower x and Q 2 in which resummation will be 
needed. 
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